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We critically review the role of cosmological moduli in determining the post-inflationary history 
of the universe. Moduli are ubiquitous in string and M-theory constructions of beyond the Standard 
Model physics, where they parametrize the geometry of the compactification manifold. For those 
with masses determined by supersymmetry breaking this leads to their eventual decay slightly before 
Big Bang Nucleosynthesis (without spoiling its predictions). This results in a matter dominated 
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observations of the Cosmic Microwave Background and the growth of large-scale structure. Given 
progress within fundamental theory, and guidance from dark matter and collider experiments, non- 
thermal histories have emerged as a robust and theoretically well-motivated alternative to a strictly 
thermal one. We review this approach to the early universe and discuss both the theoretical chal¬ 
lenges and the observational implications. 
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I. INTRODUCTION 

The success of Big Bang cosmology in predicting the correct abundances of the light elements and 
the existence of the Cosmic Microwave Background (CMB) rests upon the assumption that the universe 
began in a hot, dense state and as the universe expanded it cooled. This thermal history for the early 
universe can then be used to try to understand the primordial origin of dark matter particles, baryons, 
and the matter / anti-matter asymmetry. This picture is simple, elegant, and robust, but it can not be 
complete. Indeed, the primordial density perturbations - indirectly detected through the temperature 
anisotropies of the CMB and required for the growth of the large-scale structure of the universe - can 
not be causally generated in a strictly thermal universe. Instead, inflation provides the most compelling 
model for their origin, and during this epoch the universe was not in thermal equilibrium. Moreover, 
following inflation the universe must reheat providing the hot big bang. This process can occur far from 
equilibrium and result in an effectively matter dominated epoch that could last up to the time of Big 
Bang Nucleosynthesis (BBN). 

Given the importance of realizing inflation and reheating, model building has had a rich and long 
history (for recent reviews see [TJ 12])- ‘Chaotic’ inflation provides one of the simplest examples, where 
inflation is of the so-called large field typ^jwith the inflaton evolving over super-Planckian distances in 
the field space. Some of these models can lead to detectable gravitational wave signals, but they are also 
sensitive to Planck scale physics and dangerous quantum gravity corrections. As a result, a large effort has 
been put forth to understand how to construct models in the controlled framework of string theory. One 
lesson from these studies is that typically additional symmetries (such as Supersymmetry (SUSY)) and 
additional scalar fields (moduli) are generic and necessary ingredients for constructing realistic models. 
The former are necessary for protecting the inflaton potential from dangerous UV correction^ whereas 
moduli are a ubiquitous prediction resulting from both the presence of these new symmetries, as well as 
the string construction itself. If unstable these moduli can spoil inflation, or lead to undesirable late¬ 
time effects such as large entropy production after BBN or fifth force violations of General Relativity. 
Thus, although string theory provides a framework to construct realistic models of inflation, it appears 
to introduce a new nuisance - the Cosmological Moduli ProbleuJ^] @HS| ■ 

Significant progress has been made on moduli stabilization both during and after inflation (for reviews 
see muni). Many of the moduli have been shown to be stabilized with string scale masses by accounting 
for a combination of effects resulting from the presence of branes, fluxes, and/or strong hidden sector 
dynamics - all of which are expected ingredients of the theory. However, in explicitly constructed examples 
it is also a common result that many of the moduli will remain parametrically lighter than the string 
scale. The resulting cosmology differs from that of a strictly thermal post-inflationary history. These 
fields are displaced from their low-energy minima in the early universe and undergo coherent oscillations, 
mimicking a matter dominated epoctj^prior to BBN. The fields typically decay late through gravitational 
strength couplings, and the universe reheats (again) via the production of relativistic Standard Model 
and other light (beyond the Standard Model) particles. Depending on the masses and couplings to other 
fields, such as dark matter and baryons, this can lead to a new expectation for the post-inflationary 
cosmology prior to BBN. 

The rest of the paper is structured as follows. In Section [II] we review the cosmological moduli problem 
and briefly comment on its appearance in effective supergravity descriptions of string / M-theory. In 


1 Small field models of inflation typically require additional fine-tunings of initial conditions to enter the attractor regime 
(see e.g. [?i] and references within.). This tuning is in addition to that required for ensuring the flatness of the inflationary 
potential. 

2 In string theory, large field models with approximate shift symmetries providing the flatness of the potential can be 
constructed, but additional symmetries - such as SUSY - are required to protect this potential from further corrections. 

3 The cosmological moduli problem made its first appearance in the form of the Polonyi problem in the earliest versions 
of spontaneously broken supergravity. This was the observation that hidden sector models of supersymmetry breaking 
with gravity mediation often contain scalar fields with weak scale masses and gravitational couplings. Just after inflation, 
these fields would behave like non-relativisitic matter and continue to dominate the energy density of the universe until 
after BBN. 

4 This is strictly true only if the mass term in the potential gives the dominant contribution, otherwise the cosmological 
scaling of the energy density is determined by the dominant term in the potential nn. 
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Section III we discuss the effect of cosmological moduli on the thermal history of dark matter. We look 
at the non-thermal WIMP miracle and its phenomenological status, as well as axions from string theory 
and their status as dark matter candidates within non-thermal histories. In Section m we study the 
impact of cosmoligical moduli on theories of the baryon asymmetry of the universe, critically reviewing 
the Affleck-Dine mechanism and the baryon-dark matter cosmic coincidence. In Section ^ we study 
cosmological moduli in the context of the CMB and Large Scale Structure, looking at implications for 
the matter power spectrum, effective number of neutrino species, and isocurvature constraints. We end 
with a section on the challenges and future prospects of these scenarios. 


II. STRING THEORY, MODULI, AND A NON-THERMAL UNIVERSE 


A. Moduli Decay and Another Big Bang 


To see how a non-thermal history arises from moduli, we consider a single modulus^] which enjoys a 
classical shift symmetry a —y a + c. At low energy we expect the symmetry to be broken, resulting 
in a mass for the field m a . Moreover, if the presence of this new symmetry is tied to addressing the 
electroweak hierarchy problem (e.g. in SUSY or extra dimensional approaches) then we expect m<j <C m p 
with m p = 2.4 x 10 18 GeV the reduced Planck mass. The non-thermal history arises from the observation 
that there is no a-priori reason why the modulus er should initially begin in its low-energy minimum. 
Both thermal and quantum fluctuations can displace the scalar during the early evolution leading to a 
matter dominated phase PS¬ 
As an explicit example, consider string moduli in SUSY based approaches where on general grounds 
we expect that the shift symmetry of the modulus should be broken by both the finite energy density of 
inflation (which spontaneously breaks SUSY at the high energy) and quantum gravity effects [13] , These 
considerations imply additional contributions to the effective potential in the form of a Hubble scale mass 
and a tower of non-renormalizable operators, 


A Vi = -ci if I <7 2 + 


Cn 4 +n 

M n + 2 


+ . . . , 


(1) 


where in the absence of special symmetries c\ and c n are expected to be positiv^] order one constants, 
Hi is the Hubble rate during inflation, and M is the scale of new physics, e.g. quantum gravity. As this 
will be important later, we note that the dimension of the leading irrelevant operator that lifts the flat 
direction is model dependent as well as the scale of new physic^] Because the contributions ([l]| are the 
dominant terms in the potential during high scale inflation (Hi > m a ), this implies the minimum of the 
field at that time will be 


(°) 



( 2 ) 


Whereas, later the low-energy SUSY breaking mass will give the dominant contributions to the potential, 
and this minimum is near (a) ~ 0. This initial displacement from the low energy minimum eventually 


5 The case of non-thermal histories with multiple moduli was considered in m- Generically, string / M-theory theories 
have many moduli with various masses, and the lightest is typically the relevant one. 

6 One can argue rather generally for the positivity of c n based on causality, unitarity, and demanding that the model have 
a UV completion (although a few known counter-examples exist in gravitational systems). We refer the reader to na 
for details. On the other hand, the choice of c\ positive (leading to a tachyonic mass) is chosen to realize the general 
expectation that the high energy and low energy minima of the field would not generically be expected to coincide in a 
gravitational and/or time dependent background m- Ideally one would like to calculate this in an explicit, UV complete 
model. 

7 This is a generic expectation in string theories where new thresholds before the Planck scale are common place. Exam¬ 
ples include both the compactification (Kaluza-Klein) scale Mkk and string scale M s where M <C M p i is required for 
consistency of the effective theory m■ 
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leads to coherent oscillations of the field and the formation of a scalar condensate. The amplitude of 
the oscillations is determined by the initial displacement cr* ~ (<j), where for the potential ([I]) we have 
< 7 * ~ M (Hi/M) 1 ^ n+1 \ Hubble friction ceases and oscillations set in when the expansion rate satisfies 
H ss m a <C Hj. Because the oscillations scale like matter, they dilute more slowly than the primordial 
radiation (produced during inflationary reheating). Depending on the initial value cr*, the energy stored 
in the moduli may quickly come to dominate the energy density of the universe. At the time oscillations 
begin f osc ss H~ 3 ss m” 1 the initial abundance is given by 

P*(to sc) = \rnlal, (3) 

and once the oscillations become coherent (which typically takes less than a Hubble time) they will scale 
as pressure-less matter m with p a ~ m 2 cr 2 /a(i) 3 . The universe remains matter dominated until the 
field decays. Because the field is a modulus we expect it typically to be gravitationally coupled to other 
particles and so its decay rate is 


r a 


c m 3 
27r A 2 ’ 


(4) 


where we expect A ~ M p i and c depends on the precise coupling in the fundamental Lagrangian, but 
typically takes values of at most 0(100). Most of the field decayfj at the time fdecay ~ H ~ 1 ~ T” 1 
and we expect it to decay democratically to Standard Model particles and their super-partners. Any 
heavy super-partners produced will typically decay rapidly into the lightest SUSY partner (LSP), which 
is stable and could provide a viable dark matter candidate (more on this in Section III). 

In addition to dark matter, light Standard Model particles that are produced will thermalize and ‘re¬ 
heat’ the universe for a second time (with inflationary reheating occurring early and at high temperature). 
The corresponding reheat temperature is given by T r ~ g * 1 ^ y jV (T M p i or 


T r = c 1 / 2 



\ 3/2 rp 

50 TeV / BBN ’ 


(5) 


where Tbbn — 3 MeV and g„ is the number of relativistic degrees of freedom at T r , and this temperature 
must be larger than around 3 MeV to be in agreement with BBN light element abundances (18]. 

A comparison between the standard thermal history and that of the non-thermal history appears in 
Figure [T] for an interesting case where the modulus decays slightly before BBN and inflationary reheating 
occurs at a high temperature. We now turn to the question of what moduli masses (and so reheat 
temperatures) are naturally expected from string / M-theory approaches. 


B. Fundamental Theory and Moduli Masses 

For an interesting departure from a thermal post-inflationary universe the mass of the lightest modulus 
should be in a window between 0(10) — 0(1000) TeV. For masses much further above the upper bound 
little departure from a thermal universe is expected, whereas the lower bound is required for consistency 
with BBN as seen from Q. As a phenomenological approach non-thermal cosmologies have been consid¬ 
ered for their implications for dark matter (see [T9] for a review). However, lacking theoretical motivation 
these models are arguably somewhat exotic. 


In most of the literature the moduli decay is treated as instantaneous at tdecay- However, this approximation can be 
misleading. For example, in the case where the radiation energy density significantly drops below the moduli energy 
density (p a p r ), the continuous (even though small) amounts of particle decay before reheating can lead to changes 
in the scale factor - temperature relation because the entropy is not conserved. We refer the reader to m for further 
discussion. 
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FIG. 1: The left timeline represents a thermal history for the early universe where dark matter WIMPs 
are populated in the thermal bath that emerges shortly after after inflationary reheating. The right 
timeline represents a moduli dominated, non-thermal history, where the universe evolves as a matter 
dominated universe up until the decay time which must be before BBN. 


The key question is - should one generally expect a modulus in this mass range? In string / M-theory 
approaches, inextricably connected with the question of moduli masses is the question of the gravitino 
mass. This is because moduli stabilization plays a central role in the nature of supersymmetry breaking. 
The potential of the moduli determines their F-terms dynamically, and along with the Kahler metric of 
the matter sector, determines the strength of gravity mediation. To extract the mass of the modulus, it is 
sufficient to write down its scalar potential in the supergravity limit (given in terms of its superpotential 
and Kahler potential, W(a) and K(a)) 


V = 


e K/m 


2 

P 



3|W| 2 \ 
m 2 J 


( 6 ) 


One requires, at the extremum, that d a V = 0 and d a<J V > 0. Moreover, V ~ 0 is required to obtain an 
almost vanishing cosmological constant. We note that the gravitino mass is defined as to 2 / 2 = (e K \W\ 2 ) 
where we work in units of m p = 1 for simplicity. 

This set of conditions is restrictive enough to give a relation between m a and m 3 / 2 , without going 
into the details of visible sector model-building. Since the modulus mass matrix is positive definite by 
assumption and both terms in (Jg]) are comparable to m 2 / 2 due to the condition V ~ 0, one obtains that 
the smallest eigenvalue of the mass matrix should be 0 ( 013 / 2 )- For a single modulus, this trivially implies 

However, the modulus must be redefined to the physical basis with canonical kinetic term. An example 
Kahler potential is given by K(a) = — 21nV, where V is the volume form whose expression in terms of 
the moduli depends on the compactification manifold through the intersection numbers. In a basis with 
canonical kinetic terms, the lightest modulus will generally have a non-trivial profile along different field 
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directions. However, for many settings that provide realistic phenomenology, the lightest modulus is also 
the largest, and one can approximate V ~ er a , where a ~ 0(1). This is the case in single modulus 
examples like KKLT [20], as well as the cases of Large Volume Compactifications in Type IIB m and 
G 2 compactifications of M-theory [22]. In that case, the canonical normalization yields m a ~ ( a)d ao V , 
where by abuse of notation a is now the physical field. One thus obtains 


m a ~ (cr)m 3 / 2 . 


(7) 


Usually, one finds (a) ~ 0(10) — 0(100) in the supergravity regime, 
stabilized supersymmetricalh 9 
circumstances one expects 


so that DW ~ 0, then (a) 


log ( 


m 3/2 


For example, if the modulus is 
~ 40. Then, under quite general 


m<j ~ 0(1 - 100 ) x TO 3/2 . ( 8 ) 

For TO 3/2 ~ 10 — 100 TeV, as one would expect in models that at least partially address the electroweak 
hierarchy problem, we see that m a falls within the mass range where a non-thermal history results. 

We would like to emphasize the interesting result that the lower bound on the mass required to 
address the cosmological moduli problem provides further motivation for a little hierarchy 
between the electroweak scale and the scale of new scalar physics. Indeed, the discovery of a 125 
GeV Higgs at the Large Hadron Collider (LHC) - and nothing else - appears to have left the relevance 
of supersymmetry (SUSY) in question, particularly as a mechanism for stabilizing the hierarchy between 
the Electroweak and Planck scales. Standard Model Light superpartners (‘Natural’) SUSY still remains 
possible, but with some cost of increasing the level of complexity of models (see e.g. |21|). Within the 
Minimal Supersymmetric Standard Model (MSSM), there are a number of alternative ways to reconcile 
SUSY with the data, including models of Split-SUSY [25ff27| . or simply by accepting that some fine- 
tuning may just be an accident of naturtf 9 ] Regardless of one’s viewpoint, it seems that if low-energy 
SUSY will prevail it will require the existence of a new scale at around 10 — 100 TeV. 

It is then remarkable that in these string / M-theory based approaches with moduli that the scale 
required to avoid the cosmological moduli problem is of the correct order of magnitude for realistic 
particle phenomenologwith a 125 GeV Higgs. 

Another interesting result of this string-based framework is that the reheat temperature is not a 
free parameter, but a consequence of the hierarchy between the electroweak and Planck 
scale (determined by h? susy = m 3 / 2 mp), which also helps determine the SUSY breaking masses of other 
sparticles in the theory. In both supergravity and string motivated approaches, the key lesson is that 
the reheat temperature is intimately connected to other aspects of the theory and not a free and tunable 
parameter. Given that gravitationally coupled scalars or moduli are generic in high energy UV comple¬ 
tions of the Standard Model and in no sense exotic, we see that non-thermal histories are a feasible and 
robust possibility, and perhaps an inevitable one. Given this motivation from theory, we now turn to the 
observational consequences of non-thermal histories. 


III. DARK MATTER GENESIS 

The origin of dark matter as a thermal relic is a simple and predictive scenario for its creation, which is 
also independent of the prior thermal history of the Universe. However, as we have argued so far, one is 


9 The relation between moduli and gravitino mass in this case was first noted in m- 

10 It is noteworthy that even defining the level of fine-tuning can be an issue, see e.g. m • 

11 We note that already in 2006 this mass range for the lightest scalars and gravitino was noted in m long before the Higgs 
discovery and LHC bounds favoring heavy scalar superpartners. 






led on compelling theoretical grounds to a re-evaluation of the history of the Universe prior to BBN. This 
necessitates looking afresh at the consequences of non-thermal histories on dark matter phenomenology, 
a topic to which we turn in this section. In our choice of dark matter candidates, we will discuss both 
Weakly Interacting Massive Particles (WIMPs) and axions. 


A. WIMPs 

If the primordial universe is thermalized, massive long-lived particles may freeze-out with a final abun¬ 
dance determined primarily by their mass and annihilation cross-section [5(1] . This is the basis of thermal 
WIMPs, which assumes a weak-scale cross-section, a and ( av)th — 10 -26 cm 3 /s, where v is the typical 
velocity. Alternatively, non-thermal dark matter is produced via the decay of heavier particles into a 
long-lived final state and does not require thermal equilibrium [121I5TH55] (for reviews, see [191 1571155] , 
and for recent related work |3SfiU). Dark matter models are constrained by both direct detection exper¬ 
iments and searches for astrophysical signals generated by their annihilation products. Non-thermal dark 
matter can have a higher annihilation cross-section than thermal dark matter so astrophysical signals are 
potentially stronger for these scenarios, particularly in indirect experiments such as FERMI and AMS-2 
0M3- We now review the thermal and non-thermal WIMP paradigms. 


1. Thermal WIMPs 

The standard picture of the cosmic history assumes that the universe is radiation dominated prior to 
Big Bang Nucleosynthesis (BBN) and that dark matter is created from the thermal bath created at the 
end of inflation. The evolution of dark matter is described by the Boltzmann equation 

h x = -3 Hn x - (crv) [n 2 x - n 2 eq \ , (9) 

where (av) is the thermally averaged cross-section, nx is the number density, and n eq is the number 
density of the species in chemical equilibrium, i.e. XX o 77, where 7 is a relativistic particle. In the 
early universe with T > mx, the dark matter is relativistic and in thermal equilibrium with the plasma, 
and its density is determined by the full solution to ([ 9 ]). When T < mx, it becomes Boltzmann suppressed, 
until finally the expansion rate becomes larger than the annihilation rate, leading to thermal ‘freeze-out’. 
The dark matter density at this time is sensitive only to its mass and annihilation cross-section. 

Thus, the thermal paradigm is simple, predictive, and more or less independent of the underlying 
microscopic physics. Moreover, there is the suggestive fact that the observed relic density is satisfied for 
(a x v) « 10 _26 cm 3 -s _1 (or cr x « 1 picobarn), which corresponds to weak scale interactions. The weak 
scale is precisely where new physics is expected to appear, on very different physical grounds. Plausibly, 
then, dark matter is a Weakly Interacting Massive Particle (WIMP) that is there as a part of new physics 
at the weak scale, motivated by issues like the gauge hierarchy problem. The thermal WIMP miracle is 
precisely this combination of simplicity and serendipity. Since the weak scale is also the energy scale that 
we are currently able to probe experimentally, the thermal WIMP is by far the most actively studied 
dark matter candidate. 


2. The Non-thermal WIMP Miracle 

Despite its appeal, we saw in Section [TT| that there are compelling reasons to expect that the history of 
the universe prior to BBN was not as simple as the thermal WIMP paradigm requires. It is worthwhile 
to revisit the assumptions underlying the arguments above. They are: (i) WIMPs reached chemical 
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equilibrium in the early universe; (ii) the universe was radiation dominated at the time of freeze-out; 
(Hi) following freeze-out, there was no dilution of the relic density due to entropy production from any 
sources, and conversely, there was no further production of dark matter from late-time sources. 

The above assumptions look difficult to realize from a top-down perspective (and are in increasing 
tension with data from the bottom-up perspective, as we will discuss in Section III A 4 ). Indeed, as we 
discussed in Section ^assumptions such as WIMPs achieving thermal equilibrium, the absence of entropy 
production following thermal freeze-out, and strictly a radiation dominated universe prior to BBN seem 
to be in conflict with top-down approaches to model building. In particular, the presence of moduli imply 
violations of all these conditions. 

This by itself does not mean, however, that the thermal WIMP paradigm itself is disfavored, although 
a purely thermal history of the universe prior to BBN is. Moduli may decay early enough (above the 
freeze-out temperature of WIMPs) to render a prior non-thermal history irrelevant as far as low-energy 
observables are concerned. The entropy dilution from such very early decay is likely to erase the prior 
dark matter density and produce it afresh (by reheating the modulus into the visible sector generally). 
The dark matter thus produced will then have time to thermalize and give the usual thermal WIMP 
story. 

A non-thermal WIMP is obtained if a modulus decays after the freeze-out temperature but before 
the onset of BBN. This translates into a window in terms of the decay width, and hence in terms of 
parameters of the fundamental theory (such as the modulus mass and overall decay constant). Leading 
up to ([8]) we have argued that it is natural for moduli to fall within this mass window, and in fact the 
explicitly studied compactifications in the literature have such moduli. 


3. Reheating Again and Dark Matter Production 


Moduli decay dilutes any previous population of dark matter by a factor f l c dm —► £l c dm(T r /Tf) 3 , which 
is typically large enough to render it irrelevant. This scaling can be understood from the fact that the 
temperature goe^jas T ~ l/a(t), while volumes scale as ~ a(t) 3 (where aft) is the scale factor). 


The moduli decay is expected to produce dark matter and as a result there are two options: 


(1) If the WIMP number density exceeds the critical value n% = 


H 


(cr x v) 


then the WIMPs will 


T—T r 


quickly annihilate down to this value, which acts as an attractor \TZ\ . The fixed point value is evalu¬ 
ated at the reheating temperature instead of the freeze-out temperature. This results in a parametric 
enhancement of the relic density 


VZ = Zrn ■ ( 10 ) 

Since the relic density scales inversely with the annihilation cross section, this implies that WIMP can¬ 
didates with annihilation cross sections larger than the canonical value by a factor (Tf/T r ) can give the 
correct relic density in a non-thermal history with the above compensating factor. Examples of such 
WIMPs are the Winos and Higgsinos of the minimal supersymmetric extension of the Standard Model 
(MSSM). For future reference, we will call this option the ’annihilation’ scenario. 


If there is significant decay of the modulus during the moduli dominated epoch entropy is not conserved and this alters 
the temperature / scale factor relation HZI- In this review we will work in the instantaneous decay approximation and 
refer the reader to 1481 for a more rigorous treatment. We note that for our considerations here we do not expect any 
significant qualitative differences. 
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In more detail, the standard expression for the thermal relic density given by 

m x \ 


V T dm h 2 = 



27t\/T0 \p c m. 


9* (vv)T f/ 


( 11 ) 


Using (10 1 , (111 and (12 1 one can estimate the relic density in non-thermal dark matter as 


* 0.76 


~ 0.10 



m x 


gV 2 (av)m v T r / 


/ m x \ /10.75\ 1,IJ /3 x 10 24 cm 3 /s 
V100 GeV / l g* ) V 


(av) 


100 TeV \ 3/2 


m a 


( 12 ) 


where in the last line we used the entropy density today so = 2.78 x 10 8 p c /h 2 with h the Hubble parameter 
in units of 100 km/s/Mpc and we chose some fiducial values with < 7 * = 10.75 the number of relativistic 


degrees of freedom at the time of reheating. We note that unlike the thermal result (12) depends on both 
the properties of the dark matter (mass and annihilation rate) and on the mass of the decaying modulus 
(m CT ). As discussed in Section [II] the scalar mass (and so reheat temperature) is not a free parameter, 
but related to the gravitino mass. The mass of the scalar (and so the relic density of dark matter) 
is controlled by the need for SUSY to generate a hierarchy between the electroweak and Planck scale. 
We have chosen a fiducial value for the annihilation rate that yields roughly the right amount of dark 
matter for the hierarchy set by the choice of low-scale SUSY breaking. The cross-section is three orders 
of magnitude higher than expected with a thermal history with important experimental consequences, as 
discussed next in Section Hll A 41 

(2) The second option we will call the ‘branching scenario’. In this case, the WIMP number density 
produced from moduli decay is below the fixed point value. This implies that there is no further WIMP 
annihilations and the final number density is that produced from the decay, i.e., n x ~ B x n<j, where B x 
is the branching ratio for scalar decay to WIMPs and n a is the number density of the scalar condensate. 
With an appropriate choice of n x (that very much depends on the microphysics), this option can accom¬ 
modate WIMPs whose annihilation cross section is higher and (crucially) also lower than the canonical 
value. Note that the latter case (a pure Bino of the MSSM) cannot be accommodated in any thermal 
setting without substantial fine-tuning of the spectrum. 


4. Phenomenological Aspects and Future Probes 

We have seen that in theories with moduli and TeV scale SUSY breaking that a non-thermal history 
is attractive. When considering the associated phenomenology, it is non-trivial to avoid this conclusion 
without introducing fine-tuning. For example, one could raise both the gravitino and modulus mass, 
re-introducing the gauge hierarchy fine-tuning, or one could decouple the modulus-gravitino correlation 
in |7]) making the modulus much heavier as in racetrack models [19l l5l?| . This introduces large fine-tuning 
that is not associated with the cosmological constant. 

From a bottom-up perspective, one of the most striking aspects of the thermal WIMP paradigm is how 
difficult it is to realize even with the large number of parameters within the MSSM. Taking Neutralinos as 
the WIMP candidates, the Wino and Higgsino - which annihilate very efficiently - must be very heavy (2.5 
TeV and 1 TeV, respectively) to give the correct relic density through thermal freeze-out [53 . Whereas a 
thermal Bino typically gives a relic density that is too high, unless its annihilation cross section is boosted 
through major co-annihilation effects [51] or mixing with the Higgsino [533. Both of these options are 
fine-tuned and the latter is highly constrained by direct detection [55;, 156]. However, within non-thermal 
histories Wino or Higgsino dark matter with masses of 0(100) GeV can be accommodated within the 
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annihilation scenario discussed above. Moreover, Bino dark matter with a non-thermal history following 
the ’branching scenario’ can satisfy the relic density constraint, without relying on any other low-energy 
ingredient. Thus, an experimental confirmation of a pure gauge eigenstate of the MSSM as the dark 
matter of the universe would seem to suggest that the universe underwent a non-thermal history, since 
otherwise one would need to accept a large amount of fine-tuning. Direct detection limits are currently 
unable to probe such pure eigenstates, whose scattering cross-section off of Standard Model nuclei is 
typically ~ 0(1O -47 ) cm 2 , although future facilities like XENON1T and LZ will reach the required 
sensitivity up to WIMP masses of ~ 1 TeV Ell- 

Indirect detection provides an enormous amount of information by putting upper bounds on the dark 
matter annihilation cross section. A joint analysis of seven Milky Way dwarf galaxies using a frequentist 
Neyman construction and Pass 7 data from the Fermi Gamma-ray Space Telescope undertaken in [581 [5H 
excludes generic thermal WIMP candidates with mass less than ~ 40 GeV annihilating to bb final stated 
Current limits using dwarf spheroidals data and the PASS-8 software from Fermi-LAT collaboration rule 
out the canonical WIMP for masses below ~ 100 GeV. WIMPs in this mass range thus over-satisfy the 
relic density in a purely thermal history. Since multi-component dark matter is not an option in this case, 
the only way WIMPs lighter than ~ 100 GeV can constitute the dark matter of the universe is if they 
followed a non-thermal ’branching scenario’ cosmological history [60) . As this exclusion reach becomes 
larger, it forces more and more potential candidates squarely into the non-thermal framework. 

On the other hand, candidates like Winos and Higgsinos with large annihilation cross section are now 
severely bounded by a combination of continuum and line searches from FERMI data mm- Continuum 
photons arise from the tree level annihilation process X°X° —► W + W~, where x° denotes the WIMP. 
Dwarf galaxy data rules out Wino dark matter masses up to around 400 GeV, while galactic center data 
rules it out up to around 700 GeV for either NFW or Einasto profiles [63]. These bounds become stronger 
(weaker) if one considers steeper (core) profiles. Recent discoveries of optical signatures of black holes in 
dwarf galaxies indicates that low-mass galaxies can host intermediate massive black holes. Accounting 
for this without spoiling the measured velocity dispersion of stars, the constraints on the dark matter 
annihilation cross section may become orders of magnitude stronger m , ruling out thermal Wino dark 
matter masses up to 1 TeV. Similar exclusion bounds hold for any general scalar dark matter candidate 
that is a triplet under SU( 2) [65] . 

CMB measurements are another source of constraint for non-thermal histories. We reserve a full 
discussion of this for Section [V] here we briefly mention how CMB observations can lead to constraints 
on the self-annihilation rate of non-thermal dark matter. These constraints arise because energy injection 
from increased dark matter annihilation can alter the recombination historyp 4 ] leading to changes in the 
temperature and polarization power spectra of the CMB [66] . We have seen that non-thermal dark 
matter can have a larger annihilation rate than thermal WIMPs, and so this would lead to stronger CMB 
constraints from the additional energy injection. The 2015 results from the PLANCK mission m imply 
a constraint on non-thermal WIMPs with masses ~ 100 GeV of f e ff(crv ) < 4 x 10~ 25 cm 3 s _1 . The 
constant / e // is a model dependent parameter that depends on the mass of the dark matter and captures 
the efficiency of the annihilations - it typically takes values / e // ~ 0.1 - 1. Although this parameter does 
introduce model dependence and some uncertainty into the constraints, with the 2015 Planck release it 
is already possible to place reasonable constraints on the annihilation rates of low-mass thermal WIMPs. 
Before reaching the cosmic variance limit, this type of constraint will be capable of ruling out non-thermal 
WIMPs completely when they have electroweak scale masses and comprise all of the cosmological dark 
matter. 

Current collider bounds on Winos and Higgsinos are rather weak, given their low production cross- 
section at a hadron collider. Approaches to study the dark matter parameter space include effective 


13 The exclusion bounds are similar for other final states. 

14 This is the period during which the universe cools to a level where CMB photons can free-stream and no longer interact 
with electrons on average. 
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operators [55] and simplified models m, with a mono-X observational signature, where X may denote 
mono jet m, mono-photon eh, mono -Z mm , mono -W m , mono-Higgs, mono-6 Ea, or mono-top 
E3- a particularly promising approach to the study of dark matter is the use of vector boson fusion 
processes [76H52] . The current exclusion prospects for pure Winos and Higgsinos are the following [53] : 
taking a background systematic uncertainty of 1%(5%), with 3000 fb” 1 of data, the 14 TeV LHC is 
sensitive to Winos of 240 GeV (125 GeV) and Higgsinos of 125 GeV (55 GeV). A future 100 TeV collider 
with 3000 fb” 1 of data would exclude Winos of mass 1.1 TeV (750 GeV) and Higgsinos of mass 530 GeV 
(180 GeV) taking 1%(5%) background uncertainty. 

The above limits are in the case when everything except the lightest neutralino dark matter is decoupled 
from the low-energy spectrum. In specific models of string/M-theory compactifications, other light fields 
may be present. For example, in cases where other species such as a neutralino NLSP or gluino are 
accessible, the dark matter exclusion limits would depend on the final states and the decay modes of 
the NLSP or the colored state to the dark matter candidate. This is clearly more model-dependent. 
We briefly describe the low-energy spectrum and collider prospects of three of the most well-studied 
compactification schemes with non-thermal histories: mirage mediation models based on KKLT, Large 
Volume Scenarios, and the G2-MSSM. 

In the case of mirage mediation models based on KKLT compactification, the dark matter candidate 
may be the Bino, Wino, or Higgsino, depending on the ratio (commonly denoted in the literature by a) 
of the gravitino mass and the universal gaugino masses at the GUT scale. The mass spacing between 
the gauginos is determined by the same quantity, which also gives the relative strengths of the moduli 
and anomaly-mediated supersymmetric contributions. The rest of the spectrum depends on the modular 
weights of the visible sector. The collider and dark matter phenomenology of these models has been 
studied most recently in the post-Higgs era in [SJ and previousljj^] in [571 [HE]- The canonical case of 
a = 1 results in a spectrum with decoupled Bino dark matter, which is largely inaccessible in colliders 
and may be probed in future direct detection experiments. However, depending on the model parameters 
(modular weights) a gluino with mass ~ 2300 GeV may also be obtained in the spectrum in which case 
a discovery may be made with 300 fb -1 of data at LHC-14. Changing the uplifting sector leads to other 
values of a which result in either decoupled Higgsinos or Winos with collider prospects described above. 

In the case of Large Volume Scenarios, a wide range of phenomenological possibilities are present 
depending on the local construction of the visible sector (specifically the matter Kahler metric) and the 
way to achieve a de Sitter vacuum, as found in [8jD]. In particular, two classes of spectra were obtained 
by these authors. The first is a split-SUSY spectrum with a decoupled scalar sector, whose collider 
prospects have been discussed earlier. The second is a CMSSM/mSUGRA scenario, which has been 
studied in detail in m- There it was shown that the confluence of collider bounds (LEP, LHC), CMB 
constraints, direct detection (LUX, XENON100, CDMS, IceCube) and indirect detection (Fermi) bounds 
result in an allowed parameter space that is characterized by decoupled Higgsino-like dark matter with 
a mass around 300 GeV. 

In the case of the G2-MSSM, the collider prospects have been studied most recently by the authors of 
eh. For a benchmark spectrum with a gravitino mass of m 3 / 2 = 35 TeV, the gluino mass is expected 
to be about 1.5 TeV, while the Wino (Bino)-like gaugino mass is about 614 (450) GeV. Three and only 
three production channels should be discoverable at LHC-14: pp — > gg, pp —► X2X1 and PP \i Xi 
where \i and X2 are respectively Bino and Wino-like. The expected signature of the X1X1 channel is 
X + X~ —■*■ W + W~ + MET. The X2X1 production channel gives the final state x^xf ~t ^ + MET, 

which should be quite a clear channel at the LHC. Given the dark matter mass of 450 GeV, a 1.5 TeV 
gluino is expected to be discoverable at the 5cr level at LHC-14 given 300 fb -1 of data. 


15 These models have also been studied in the context of S'O(IO) Yukawa-unified models l85i [861 . 



13 


B. Axions 


In addition to WIMPs, axions provide another well-motivated particle candidate for the cosmological 
dark matter. The QCD axion provides a possible solution to the Strong CP problem [92I1M] . whereas 
string axions arise ubiquitously in string / M-theory constructions upon compactification to four di¬ 
mensions [SSHSiEj. In both situations, axions arise as pseudo-Nambu-Goldstone Bosons (PNGB) of a 
spontaneously broken U{ 1) Peccei-Quinn (PQ) symmetry. The observational consequences of the axion 
depend primarily on its mass m a , the axion decay constant f a , and the initial misalignment angle 9 a . 

For the QCD axion, the mass is related to the decay constant as the axion enjoys a shift symmetry to 
all orders in perturbation theory and its mass is assumed to result from non-perturbative QCD instanton 
corrections V(a) = A 4 CD (1 — cos (a/ f® CD )), where A 4 CD = {rn a f® CD ) 2 , and so there is a relation between 
the mass and decay constant [52H93] , 

m a ~ 6 x IQ ” 10 eV • (13) 


Lab and astrophysical constraints require /^ CD > 10 9 GeV with stronger model-dependent constraints 
possible depending on the axion coupling to photons and whether the axion is assumed to comprise all 
of the cosmological dark matter [W] . 

Whereas, for string theory axions the relation (131 need not apply as the mass can be generated by 
a number of different non-perturbative effects. Stated another way, this implies challenging modeling 
building constraints for realizing a solution to the Strong CP problem in string theory - notably QCD 
instantons need to provide the dominant contribution to the low-energy effective mass of the QCD axion. 
Yet, investigating the UV completion of the PQ symmetry seems essential given the expectation that 
global (PQ-like) symmetries are expected to be broken by quantum gravity effects |100H102] . Regardless 
of whether string theory can accommodate the QCD axion these theories do predict PNGBs that can act 
as all or part of the cosmological dark matter and may have other interesting observational implications. 


1. Axions from String Theory 


Within string / M-theory constructions, PNGB’s naturally arise from the compactification of higher 
dimensional gauge fields to four dimensions [9F, 1103] . The higher dimensional gauge invariance of these 
fields leads to an approximate shift symmetry a —> a + const, in the 4 D low-energy effective theory. This 
shift symmetry can be lifted by both quantum and non-perturbative effects that depend on the details 
of the compactification [S5j. As an example, consider axions arising from compactifications of Type IIB 
string theory. To see how axions arise in the theory, we consider the dimensional reduction of the theory 
to four dimensions by starting from the 10D action in the string frame given by m 


s 


IIB 

10 


1 

(27r) 7 Q! , 4 


d 10 xV^G 





^l 2 


(14) 


where Gmn is the ten dimensional string frame metric and H 3 = dB 2 and F 3 = dC 2 are the NS-NS 
and RR three-form fluxes, respectively, with B 2 and C 2 the corresponding gauge potentials and l/(27ra') 
is the string tension. The model independent axion Co and dilaton are combined as the axio-dilaton 
r = Cq + i/g s where g s = exp(</>o) is the string coupling and we will take Cq to be fixed and instead 
concentrate on the model dependent axions arising from the compactification of the form fields. The 
additional terms represented by dots include higher form fields such as C 4 (we will concentrate on C 2 
axions). 

The zero modes of B 2 and C 2 are independent of the co-ordinates of the compact dimensions and can 
be integrated over chosen two-cycles of the internal geometry giving rise to axions in the four dimensional 
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theory. To make this explicit, consider compactifying on a Calabi-Yau 3-fold (CT 3 ) and for the form field 
C 2 we make the ansatz [95] 


C 2 


1 

2ir 


Cl{x)ui I , 


(15) 


where the Cj(x) are only functions of the four non-compact space-time dimensions and / labels the 
two-cycle. We have introduced the basis forms oj 1 to describe the internal geometry and they obey the 
normalization condition f s u J = (2ir) 2 a'5j with the two cycles Ej giving a basis of the dual homology 
H 2 (X,h). The normalization factors of 2tt are chosen for later convenience. Making a similar ansatz for 
B 2 and using this in (fTdj) we have 


S = 


d 10 x 
(2n) 7 a ' 4 






A_ G nrn G l Pd ^bj^ J m )j 
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(16) 

where Greek indices run over the four non-compact dimensions and lower-case latin indices denote the 
compact dimensions. At the classical level, the gauge invariance of the higher dimensional gauge potential 
implies that the axions can only be derivatively coupled and so we have a shift-symmetric pseudo-scalar 
in the low energy theory. This symmetry can be broken in a number of ways, which we will discuss 
shortly. Denoting both types of axions as a/, upon dimensional reduction we find 


S 4 
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7 IJ d ll aid tl aj 


(17) 


where g is the 4 D Einstein frame metric and the 4 D reduced Planck mass is 


m 


2 

v 


2V 

(2tt ) 7 g 2 a' 


(18) 


with V = Ve/a' 3 the string fraiiHp^] volume of CY :i . 

For the RR axion the r ) 1J provide the axion decay constants and depend on the internal geometry as 

ESI 
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ij 

RR 


1 

6(27r) 9 a' 4 


J Ull A -kUJj, 


(19) 


whereas for the NS axion one gets the same result multiplied by an extra factor of gj 2 . Once the 
specifics of the internal geometry are known one can calculate the 'y IJ to find the corresponding axion 
decay constants. This is a non-trivial task, which requires a full specification of the compactification 
geometry in a Calabi-Yau manifold. In |jj5] the 7 /J were calculated in a variety of string models assuming 
compactifications sufficiently symmetric to be amenable to estimates. 


Generating a potential. Although the higher dimensional gauge invariance guarantees a low energy 
shift symmetry to all orders in perturbation theory, there are a number of non-perturbative effects that 
can generate a potential - e.g. branes wrapped on the corresponding cycle of the axion [1011105J and/or 
various instanton and condensation effects [97, lOli 1106H108j . These effects can completely lift the axion 


16 Another common convention is to instead work with the Einstein frame volume. The 10D string frame is related to the 
Einstein frame by the Weyl rescaling = exp(0/2)G|^ steln and working in units of the string length l s = 2n\/a' 

the two volumes are related by Ve = (27r) 6 V/<?s 2 . 
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from the low energy spectrum or result in a cosine-like potential just as in the case of the QCD axion, 
albeit at vastly different scales. In general, the range of axion masses expected from these effects can 
range over an entire spectrum from very heavy, to string instantons generating masses far smaller than the 
QCD axion. The large range of model building possibilities has interesting observational implications for 
many different types of experiments and the resulting theoretical picture has been termed the "axi-verse" 

m- 


Moduli Stabilization, the Saxion and a Non-thermal History. As familiar from inflationary model 
building in string theory, it is difficult to introduce a single degree of freedom without examining other 
ingredients of the fundamental theory. String / M-theory compactifications come with a number of 
additional ingredients, which is the cost for working with a UV complete approach. In particular, one 
must take into consideration the additional moduli, fluxes, and branes within a chosen framework. We 
will not review these difficulties in detail here (see e.g. [97 3 98! and references within). However, we note 
that moduli stabilization with light axions has a generic challenge as its SUSY partner - the saxion - can 
lead to a number of difficulties. For example, if the saxion is stabilized supersymmetrically (as Kahler 
moduli are in KKLT setups before the uplift to deSitter) then the saxion and axion are both expected to 
have comparable masses. Thus, if the axion is to remain light and the saxion receives a comparable mass 
then this can lead to an example of the cosmological moduli problem - discussed in Section [IT] Whether 
this is the case or if there can be a large mass splitting is somewhat model dependent. However, this must 
be addressed within a given setup to guarantee a successful axion dark matter candidate. As an example, 
for a saxion receiving a mass of ~ 50-TeV this once again alleviates the cosmological moduli problem 
and motivates a non-thermal history since the saxion would decay slightly before BBN. Depending on 
the initial conditions for the axion this may or may not lead to a successful model for axion dark matter 
and agreement with the relic dark matter abundance. 


2. Axion Dark Matter Density in a Non-thermal History 


Given motivation for axions from both string / M-theory and the Strong CP problem in QCD we now 
consider the role of the axion as dark matter in a non-thermal universe. We will not review the standard 
case of axions in a thermal universe, instead we refer the reader to existing literature j!09lilH~] . However, 
the standard calculation for the relic density of axion dark matter will hold, even in the non-thermal 
history with moduli decay, if the axion mass is such that axion oscillations begin following moduli decay 
- i.e. if m a < r CT with T a the moduli decay rate given in This typically holds true for axion masses 
m a < 10“ 15 eV, and then axion oscillations will begin following moduli decay and the critical density in 
axions today is m 


n a h 2 ~ o.i 


fa 


2.0 x 10 16 GeV 



10- 20 eV 



( 20 ) 


where fl a = p a /p c is the density of axions today normalized to the critical density and this expression is 
accurate up to 0 { 1) corrections due to effects such as anharmonicity m- We note that for the special 
case of a QCD axion (201 can be further simplified by making use of (131. As mentioned above, the result 
(201 depends on the axion mass and decay constant, as well as the misalignment angle of the axion from 
its minimum 6$. In most of the literature the initial value 9q is taken as random variable with an average 
squared value (0q) ~ 7t 2 /3. However, the actual value depends on a number of assumptions, including 
the status of the PQ symmetry during inflation |H'3I . In this review we will leave this value unspecified 
to be as general as possible, invoking the standard average when numeric answers are required. We see 
from (20) that for a decay constant in the vicinity of the GUT or string scale ( f a — 10 16 GeV) we must 
either have extremely light axions m a <C 10 -20 , or fine-tune the initial misalignment angle -C 1, in 
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order not to over-close the universe, i.e. to be consistent with cosmological observations which imply the 
abundance of axions today must be f l a h 2 < 0.1. 


Axion oscillations before moduli decay. Next we consider when the axions begin oscillating prior to moduli 
decay requiring m a > 10“ 14 eV. In this case the moduli decay and the entropy produced will dilute the 
relic abundance as first discussed in [ 1141 . For completeness, we now briefly review this calculation of the 
relic abundance as it appears in [112 . 

The equation of motion for the axion in the background of the oscillating moduli is given by 

a + 2>Ha + m 2 a = 0, (21) 

where H ~ 2/(3 1) in the effectively matter dominated epoch and m a {T ) is the axion mass which can 
now depend on the temperature ma. The energy density of the axion is p a = (a 2 + ?n 2 a 2 )/2 and once 
coherent oscillations begin (typically in much less than a Hubble time, A tjj ~ H ~ l ) then the average 
energy density evolves as (p a ) = (0 2 ) = ( m 2 {T)8 2 ) and the pressure then vanishes. We then have 


Pa = - 



( 22 ) 


where the time dependence of the mass follows from the temperature dependence and so the energy 
density evolves as 


m a (T) 

Pd PO /,\9 •> 

a(tr 


(23) 


where aft) is the scale factor (not to be confused with the axion). The temperature dependence of the 
mass implies that the energy density does not scale like matter but that the number density (n a ) does. 
Thus, if the axions start to oscillate during the moduli phase their number density compared to that of 
the moduli (which we recall also scale as ~ 1 /a(t) 3 ) is constant 


±m a (T^)f 2 a Ael ff 

m„ 


where H osc and T osc denote the Hubble parameter and temperature at the onset of axion oscillations 
t = t osc and A— f((@o)) * s the axion displacement including various finite temperature and other 
corrections which can be found as a function of the zero temperature displacement 9q m■ The amount 
of axions after moduli decay and at the time of reheating t r is 


na _ na ( Pc r \ 
S t~>tr Pa \yS(T r ) J 


(25) 


where s(T) is the entropy density and we assume that the moduli do not decay into axions (we will 
discuss other cases soon). Comoving energy is conserved during the decay so for instant reheating we can 
relate the moduli density to the radiation after reheating and p„ = p r = 7r 2 g*T 4 /30, so that 
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where we used s(T r ) = 27r 2 g s (T r )T)?/45 and g s {T r ) 
the critical density is 


g*(T r ). Using that (26) is conserved until today, 
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where Sq and p c = SHfirri 2 are the entropy and critical density today, respectively. Using p c j Sq = 


3.7 x 10 9 h 2 GeV with h the Hubble parameter in units of 100 km/s/Mpc and (24) and (261 we have 

7„ 2 A9 2 T, 

^ Cose m2 

fa V ( T, 


n a h 2 = 3.o x io 8 
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10 16 GeV 


1 MeV 
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(28) 


where £ osc = m a (T oe . c )/m a ( 0) is the ratio of the finite temperature mass to the zero temperature mass 
and is typically an order one number [ 1112 1 . This result is interesting in that as long as m a > 10” 14 
eV the resulting abundance is independent of the mass of the axion. We see from (28), for low reheat 
temperatures and axion decay constants near the GUT scale f a — 10 16 GeV, the initial displacement 
angle would require some level of fine-tuning (A0 e // -C 1) in order to be consistent with the observed 
dark matter abundance. This presents a challenge for many string and M-theory based setups containing 
axions. For example, in EH it was found that for M-theory approaches with gauge coupling unification 
that the masses of the axions would be distributed logarithmically between 10” 33 eV and about an eV, 
but that typically f a should be near the GUT scale. However, it is again important to stress that a 
knowledge of the nature of the PQ symmetry during and after inflation is important for establishing the 
likely values of the displacement [113j . 


3. Observational Constraints on Axions 

In addition to bounds on axions from the dark matter relic density, their behavior both during inflation 
and at the time of moduli decay can provide important constraints. For axions with masses comparable 
or lighter than the Hubble scale (Hi) during inflation m a < Hi axions can generate an isocurvature 
component in the primordial perturbations. Whereas, axions which are relativistic near the epoch of BBN 
can act as an additional source of dark radiation which is also constrained by both CMB observations 
and BBN. Both of these constraints can be quite stringent and depend on assumptions of the cosmic 
history. We will discuss both of these constraints and their connection with the moduli phase in Section 

EU 

Many of the constraints on axions are rather independent of whether the post-inflationary history is 
thermal or non-thermal. We will mention these briefly and refer the interested reader to the literature for 
more details - see e.g. and references within. The relevant constraints depend primarily on the mass 
range of the axions and on the strength of the coupling to photons ~ aE ■ B. In [96] . it was shown axions 
with masses m a < 10” 28 eV (and with appreciable electromagnetic coupling) lead to rotation of the CMB 
and the level of constraint will depend on the strength of the coupling - of course for axion dark matter 
candidates this coupling must be quite weak. For axions in the mass range 10” 28 eV < m a < 10” 18 
eV step-like features in the matter power spectrum can become important for structure formation on 
small scales. For larger masses in the range 10” 18 eV < m a < 10” 10 eV axions can form bound states 
with black holes leading to a spin-down of the black hole and a mass gap in the spectrum. Finally, for 
axions with masses larger than around 10” 10 eV and with a significant coupling to photons, axion-photon 
conversion can occur in strong magnetic fields resulting in signals both in galaxies and in more compact 
objects such as pulsars. 


IV. BARYOGENESIS AND THE COSMIC COINCIDENCE 

String cosmology has been most widely studied in the contexts of inflationary, axionic, and dark matter 
physics. On the other hand, the baryon asymmetry typically relies on model-building in the visible sector. 














While realistic vacua that approximate the Standard Model have been constructed in different corners of 
M-theory, it is safe to say that extensions required for baryon asymmetry have not been a model-building 
priority. This is understandable, since processes such as electroweak baryogenesis m or leptogenesis 
m are relatively independent of a particular UV completion, according to the common lore. 

There is every reason to expect that the viability of canonical processes like electroweak baryogenesis 
and leptogenesis does, in fact, have an unavoidable dependence on UV physics. The decay of moduli will 
dilute any previous population of baryons. Since electroweak baryogenesis depends on processes at scales 
~ 0(100) GeV, and moduli decay at scales ~ 0(1) GeV, this dilution factor can and should be large. 

There are then two options. The first is to rely on a high scale baryogenesis process such as Affleck-Dine 
(AD) baryogenesis which produces a large baryon asymmetry that can withstand subsequent dilution. 
The second is to implement baryogenesis at a post-sphaleron stage. We discuss these options in turn. 


A. AD Baryogenesis: UV Constraints and Challenges 


Affleck-Dine (AD) baryogenesis |T3], [15], ( ITT] relies on very general properties of supersymmetric 
theories: the existence of flat directions. An inflationary sector produces coherent oscillations along a 
supersymmetric flat direction S , and the interaction between the inflationary sector and the flat direction 
leads to dynamics in field space that ultimately breaks CP and baryon number. This interaction occurs 
through Planck-suppressed operators. The setting is effective D = 4, N = 1 supergravity, and string 
phenomenology plays a definite role. 

The relevant potential for the AD field S in the early Universe is mm- as explained below, 


V(S) = (-c H H 2 +m 2 )\S\ 2 
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where c# and a correspond to a Hubble-induced mass and A-term, whereas and A are the soft- 
SUSY breaking mass and A-term, respectively. The last term in the potential corresponds to a higher 
dimensional operator in the superpotential, W D A S n /M n ~ 3 (M is the cutoff scale where new physics 
appears). 

The curvature along S is dominated by the Hubble-induced mass term for H m so f t . If cjj > 0, the 
field sits at the origin. However, if ch < 0, the minimum is non-zero and the field tracks this minimum 
until H ~ ms- The field settles into one of the n discrete vacua given by the Hubble induced A-term. 
When H ~ ms, the new minimum develops at the origin and the field begins to oscillate around it. The 
soft A-term becomes dominant over the Hubble-induced A-term, and the field acquires an angular motion 
to settle into a new phase. The baryon number violation becomes maximal at this time and imparts an 
asymmetry to the condensate. 

The different components of (291 are non-trivial to justify from a top-down perspective. We describe 
these challenges in turn. 


(i) Negative Hubble-induced mass. The foremost requirement is a tachyonic Hubble-induced mass 
term for the field S. This requires a condition on the geometry of the moduli space. The Kahler potential 
and superpotential can be written as 

K = K{<Ti,ai) + K a p(a il o : i )S a SP + ... 

W = Wfa) + \Y afil S a ^ , 

where S denotes a generic visible sector field. The soft mass is given by 


( 30 ) 











19 


with the holomorphic bisectional curvature H[jS, a)) being given by 

H[5, cr] = - . (32) 

9aa9ss 

where R denotes the Riemann curvature tensor of the moduli space. 

Since during inflation we can take Vo TO 3 / 2 , the requirement of a tachyonic Hubble-induced mass 
reduces to 


U[S,a] < -1 . (33) 

This condition has been explored in a wide variety of contexts within type IIB string theory [TT 8 J 11T9] . 
For a non-local modulus cr, it was found that the Hubble-induced masses are generally positive, with 

ch = 2/3 , (34) 

making AD baryogenesis difficult to satisfy. If the modulus cr corresponds to deformations of the local 
cycle on which the visible sector is constructed, on the other hand, then the induced masses are also 
typically positive, although for suitable choices of K a jj(ai,ai) a tachyonic mass can be generated. If the 
visible sector construction is at a singularity, then the induced mass depends on appropriate choice of 
fluxes on the mirror type IIA manifold. Aspects of this problem have recently been discussed in |120H123] . 


(ii) Hubble-induced A-terms. To obtain successful AD bayogenesis, the presence of a Hubble-induced 
A-term is necessary. The A-term is given by 


A a pj — F l ( diY a [ 3 7 + -KjYaf 3 7 — 


(35) 


where T is the connection on moduli space, and the *-th modulus cq acquires the dominant A-term during 
inflation (we will simply call it a henceforth). 

We need the presence of a Yukawa coupling Y al 3 7 in the superpotential in (301. To obtain Hubble- 
induced A—terms during inflation, the various contributions in (|35|) should acquire a non-zero value. 


Generally, the dependence of the Yukawa coupling on the modulus sector depends on the details of 
the local model building. Consider a visible sector localized near a small four-cycle r s = ReT s . Due to 
holomorphy and the shift symmetry of the Kahler moduli, they cannot appear at any level in perturbation 
theory in W, and hence, importantly, in Y of (30). The entire dependence on Kahler moduli appears in 
the normalized Yukawa couplings, given by 


Y afh (r s ,U) = e */ 2 YaM{U) 


(K a KpK^' 


(36) 


This should only depend on local geometric data r s and complex structure moduli IA , but not the overall 
volume or other Kahler moduli. 


The first term in (35 ) requires that the Yukawa couplings depend on the modulus (cr) that acquires the 
dominant A—term during inflation. This term vanishes in general since the Kahler modulus corresponding 
to deformations of the cycle on which the visible sector is constructed is typically not the dominant source 
of supersymmetry breaking during inflation. Obtaining a non-zero value for the other two terms in (351 


depends both on the geometry as well as the specific details of inflation. A detailed study of these issues 
is merited. 

The AD baryogenesis scenario also has a number of correlated predictions coming from the presence of a 
light scalar field during the inflationary era. One important issue is the back-reaction on the inflaton. The 
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AD field (S) makes a transition from the origin to its non-zero vacuum configuration in field space when 
the Hubble-induced tachyonic mass term becomes dominant. Depending on when this transition takes 
place, there will be back-reaction on the inflationary potential, leading to features on the power spectrum, 
which will be constrained by data m- The presence of a scalar field coupling to the inflationary sector 
should also lead to constraints from non-Gaussianity. Both of these issues deserve further study. 


B. Moduli and the Cosmic Coincidence Problem 


Cosmological observations not only determine precisely the relic abundance of dark matter and baryons, 
but also imply an interesting connection between their relative amounts f ^dm/^B ~ 5 leading to what 
some have called the ‘Cosmic Coincidence Problem’. In nm it was argued that moduli decay in non- 
thermal histories might address this ‘coincidence’ - the primary feature being that both dark matter and 
baryons can be produced from the decay of the lightest modulus. Keeping in mind that moduli decay will 
dilute any initial abundance of particles, it is advantageous in non-thermal histories that AD baryogenesis 
typically over-produces baryons. Indeed in [ 125] the authors pointed out that the late decay of a Polonyi 
field could act to dilute the baryon asymmetry to an acceptable level. More recently, in [ 124] it was noted 
that because dark matter can also be produced in the decay, this creates a subtle connection between the 
two relic densities that could account for the baryon-dark matter coincidence problem. To understand 
the connection it is useful to consider the ratio of the baryon to moduli co-moving densities 



where <7o is the amplitude at the start of the moduli oscillation, and the subscript i and superscript 0 
denote the values before and after modulus decay, respectively. 

To obtain the baryon number density after decay, the dilution factor A = s a fter/sbefore must be used. 
We obtain 


K = 4^ _ (Yg\ _ 3 Tr_ /Sq\ 2 /«b\ 

A Safter Safter \ Y° ) 4 m S \CT 0 J \n S ) i 


(38) 


which demonstrates that the baryon abundance in this approach is intimately related to the ratio of the 
initial amplitudes of the AD field and the modulus. The AD fields displacement depends nontrivially on 
the dimension of the non-renormalizable operator that lifts the flat direction. Since larger n leads to larger 
displacement, and therefore larger contribution to the baryon asymmetry, one can focus on the flattest 
directions in MSSM that require the largest n to get liftecFd As showed in Ref. |126| . the flattest direction 
(one of the Q,u, e combinations) corresponds to n = 9. Assuming that the non-renormalizable operator 
is generated at the reduced Planck scale M ~ M v and taking ms ~ 10 5 GeV we find So ~ 10 16 GeV. 
For the next flattest direction (one of the d, L combinations) - which is not lifted until n = 7 - we have 
So ~ 3 x 10 15 GeV. So we can see that these directions in the MSSM naturally have amplitudes two or 
three order of magnitudes smaller than M p , and it follows that the correct baryon asymmetry can be 
obtained for T r ~ 100 MeV, ms ~ 100 TeV, and (So/oo) ~ 10~ 2 . 

The modulus will also decay and provide dark matter, and the relic density will depend on whether the 
number density produced is above (’annihilation scenario’) or below (’branching scenario’) the critical 
number density. In either case, the final baryon-dark matter relic density ratio will depend on the factors 


17 We assume all non-normalizable operators that are allowed by gauge invariance and R-parity are generated. 
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in (38) and the annihilation cross section of the dark matter or the branching ratio of a into dark matter. 
By a reasonable choice of parameters, one would obtain f /Qdm ~ 1 : 5. 

We note that similar approaches to the baryon-dark matter coincidence problem that depend on moduli 
physics have been undertaken recently |127H1321 . A more direct solution to the coincidence problem would 
result if the both the baryons and the dark matter resulted from the decay of the modulus. This idea was 
first proposed in j!33] , and it would be very interesting to see if this approach could be UV completed in 
String or M-theory. 


C. Baryogenesis at Late Times 

Another approach to baryogenesis in non-thermal histories can be implemented by a leptogenesis- 
like mechanism resulting from modulus decay. This approach invokes additional model building, but it 
implies an intriguing connection between the need for Y a = (3T r /4m CT ) ~ 10 -8 to be a small number 
and that baryogenesis requires small number densities ~ 0(1O -10 ). This small value can be driven by 
the smallness of Y a , with the remainder factor of 10~ 2 coming from branching ratios and loop factors 

[TM1H55I . 

As an explicit model, we can consider the following extension of the MSSM, containing two flavors of 
singlets N a and a single flavor of colored triplets X, X (with hypercharges +4/3, —4/3 respectively). The 
superpotential is given by 

Wextra = Aj a N a u^X + X'^d^djX H— —N a N a + MxXX . 

(39) 

The interference between the tree-level and one-loop diagrams in the N a —>• Xu c decay generates a baryon 
asymmetry, in a manner reminiscent of leptogenesis. The asymmetry per decay e is loop suppressed and 
naturally e ~ 0(0.1) for 0(1) phases and couplings. Taking Y a ~ 10~ 7 , the correct baryon asymmetry 
can be obtained for a branching of the modulus into the species N given by Br jv ~ 10 -2 . From counting 
of the degrees of freedom this branching ratio appears to be somewhat natural. 

The modulus will generally produce dark matter \ in addition to the the field N required for baryo¬ 
genesis. Considering that the number density of the dark matter produced is below the critical density 
required for further annihilation, one has for the baryon and dark matter density ratio 

Q b 1 GeV e Br^v 

o- - - x “r- ' 40 

S1 DM m x Br x 

The branching scenario requires Br x < 10 -3 . One then obtains the correct coincidence for m x ~ 0(10) 
GeV, but other values of m x can also be accommodated. 


V. GRAVITATIONAL CONSTRAINTS: CMB AND LARGE SCALE STRUCTURE 

In this section we consider possible implications of a non-thermal history on the CMB, BBN and the 
Matter Power Spectrum. Thus far, many of the observational consequences of a non-thermal history have 
relied on a connection with interactions or decays to Standard Model particles. In the case of dark matter, 
we saw that non-thermal WIMPs can have unexpected properties (e.g. a stronger self annihilation cross 
section than the thermal case) and as a result we saw that (in)direct detection of dark matter experiments 
could help constrain the reheat temperature. However, these probes rely on the assumption that dark 
matter is composed of WIMPs and have a connection to electroweak-scale interactions. 
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FIG. 2: Evolution of (normalized) dark matter contrast {8dm = 8pdm/Pdm) in a non-thermal cosmology 
for modes entering the horizon at IV ~ 5 (Left) and N ~ 7 (Right) where $o is the initial metric 
perturbation. As the mode enters the horizon it grows linearly with the scale factor and modes entering 
earlier experience more growth. After the universe becomes radiation dominated at N ~ 12, the 
amplitude of the density contrast decreases due to rapid annihilations of dark matter particles, and 
then the density grows logarithmically as expected. 


What if dark matter is not directly associated with weak-scale physics? We discussed one example 
- axions - and we saw that if these particles are light they could prove difficult to detect. If part of 
the cosmological dark matter has negligible interactions with (or decays to) the Standard Model, this 
would present a substantial challenge to identifying its particle nature. Moreover, as LHC prepares to 
probe energies in excess of 10 TeV, there is, as yet, no evidence for physics beyond the Standard Model. 
If nature has decides to play a cruel trick on us and the scale of new physics lies significantly above 
the electroweak scale, our only near-term possibility for probing new physics may rely on cosmological 
observations. In this section we present implications of the non-thermal history resulting from string / 
M-theory approaches on gravitational and cosmological based experiments. In particular, we consider 
what the early matter dominated phase could imply for the formation of primordial structures (such as 
mini-halos and primordial black holes), we also discuss implications for the matter power spectrum, and 
different ways in which the matter dominated phase can influence observations of the CMB. 


A. Implications for the Matter Power Spectrum 

1. Enhanced Structure Growth on Small Scales 

In conventional models for structure formation growth of structure only becomes significant following 
matter-radiation equality. This is because dark matter perturbations Sp/p that enter the horizon only 
grow logarithmically with the scale factor in a radiation dominated universe Sp/p ~ log a(t) ~ logf, 
whereas they grow linearly in a matter dominated universe Sp/p ~ aft) ~ f 2 / 3 . However, in non- 
thermal cosmologies we have seen that the universe is matter dominated prior to BBN. This suggests 
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FIG. 3: Evolution of the (normalized) radiation density contrast S r = Sp r /p r and velocity perturbation 
6 r for the inodes k/Hg = 0.1 (Left), k/Ho = 0.04 (Right) where <f>o is the initial metric perturbation. 
This mode crosses the horizon at N ~ 4.5 (Left), N ~ 6.5 (Right). In this non-thermal cosmology the 
universe is effectively matter dominated until N — 10.5 e-foldings after which the universe becomes 

radiation dominated. 


an early period of growth and a new scale at smaller wavelengths in the matter power spectrum. This 
possibility was first studied in the absence of annihilations in and later more generally in {45| . The 
results are summarized in Figure [2j It was found that the moduli phase would introduce a new period 
of growth for dark matter perturbations and that the size of the horizon at reheating could provide 
a new scale for determining the size of the smallest possible dark matter substructures. Indeed, for 
reheat temperatures near the limit set by BBN it was shown that a detectable level of Earth-sized, 
ultra-compact mini-halos may result. Of course, whether such structures can survive depends on a 
number of details, including whether free-streaming and kinetic decoupling effects would destroy the 
initial formation process. In [48], it was found that for typical WIMP candidates in string-based non- 
thermal cosmologies that free-streaming and kinetic decoupling would instead determine the smallest 
structures. However, as discussed in [4511136] . there are exceptions pointing to new directions for model 
building and a possible signature of the non-thermal period. One example is if the WIMPs are produced 
non-relativistically in the moduli decays (so that free-streaming is negligible) and if these particles become 
kinetically decoupled before reheating. Whether such models can be realized from fundamental theory 
remains an open question, although if such structures do survive they could lead to new expectations for 
indirect detection experiments. 

In addition to enhanced growth of dark matter perturbations, it was also found in [45, , 136] that there 
can be interesting effects on radiation perturbations. These perturbations also grow during the moduli 
dominated phase, and then following reheating those that entered the horizon will oscillated with an 
amplitude suppressed relative to the situation in a strictly thermal history. As seen in Figure [3] the 
suppression in the radiation perturbations arises at the peak of scalar decay when most of the radiation 
is created by both dark matter annihilations and decays of the moduli. The suppression of the radiation 
perturbations can lead to a new critical scale for dark matter structures if dark matter is produced 
thermally following reheating [48] 136] , since the suppression would prevent acoustic oscillations on that 
scale from wiping out structure [4571- However, again it was argued in [45] that this would require non- 
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standard WIMPs, since for the interesting case of low reheat temperatures (~ MeV) this temperature 
would typically be far below the thermal freeze-out temperature (~ GeV). 

In summary, it was found in |48L 136! that non-thermal histories do not face any challenges in regards to 
disrupting structure formation. Moreover, if we consider less typical SUSY WIMP dark matter candidates 
it way be possible in those models to distinguish between non-thermal and thermal scenarios through 
future observations. 


2. Primordial Blackhole Constraints 

Above we primarily focused on the growth of dark matter and radiation perturbations during the 
moduli phase. However, both of these components are typically sub-dominant in string-based approaches 
and instead the moduli significantly dominate the energy density. We argued above that the transfer of 
these perturbations into dark matter and radiation can have interesting effects in some cases. But one 
can ask if the growth of moduli perturbations during the early matter domination can lead to interesting 
constraints? In particular, it was realized many years ago that an early matter dominated (dust) phase 
could lead to a significant production of primordial black holes (PBHs) |f5§) . The predicted number of 
PBHs relies sensitively on the scalar tilt of the primordial power spectrum as this sets the initial amplitude 
of the perturbations and PBHs are expected to form when the relative density contrast approaches the 
non-linear regime Sp/p ~ 0(1). A matter or moduli dominated phase differs in some important ways 
compared to the typical estimate of PBH formation in a thermal history. In particular, during a matter 
dominated phase there is no pressure to prevent the collapse of sub-horizon perturbations to evolving into 
black holes. This means that unlike the thermal case, perturbations can collapse to form PBHs long after 
crossing the Hubble radius - the Jean’s radius vanishes and so there is no pressure to prevent collapse. 
However, in [ 138) (see also [139] ) it was stressed that in a matter (dust) phase the resulting number of 
PBHs will depend on the probability for perfectly spherical collapse (so that angular momentum doesn’t 
disrupt the process) and the duration of the period. In |140 , this analysis was adapted to the case of 
moduli domination (similar analysis with different motivation appeared in [ 1411 1'42| ). One important 
consideration is that scalar fields on small scales can cease to be homogeneous and this can introduce a 
gradient Dressure that Drevents collaose 11431 Il44l. However, it was shown in ED] that for the range of 
moduli masses of interest that this effect is negligible. It was also shown that the most relevant PBH 
constraint is requiring that PBH production does not exceed the critical density. Given, a priori, the 
scalar tilt of the primordial power spectrum the corresponding constraints on the reheat temperature 
(and so moduli mass) were found in [HD], There it was also shown that for a fixed duration of moduli 
domination and a specified reheat temperature (both which follow from the moduli mass), constraints on 
the primordial scalar tilt n s can be established. For high scale inflationary reheating and T r ~ Tbbn the 
constraint n s < 1.3 was found. However, given that at the time CMB modes excited the horizon (~ 55 
e-foldings before the end of inflation) the scalar tilt is constrained to be strictly red (n s < 1) one can 
question the importance of this constraint. This implies that only in models where significant changes in 
n s result would one be able to place strong constraints on the reheat temperaturq^q Examples of such 
the blue spectra on small scales have appeared in the literature, for example in models of hybrid inflation 
(see 1145 for a recent discussion of other possibilities and references). In summary, PBH constraints do 
not lead to strong constraints on non-thermal histories unless the scalar tilt becomes blue (n s < 1) on 
sub-CMB scales. 


18 This result may change given an extremely long period of matter domination, but it was found in m that even for high 
scale inflationary reheating that requiring radiation domination by the time of BBN forces the epoch to be too short in 
duration to place constraints on a red spectrum. 
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B. CMB Constraints 

CMB observations are not only sensitive to the physics present during inflation, but also that during 
the post-inflationary history. Significant work has gone into establishing predictions from inflation for 
both signatures in the temperature power spectrum and those associated with non-gaussianity (see e.g. 
[T]). In fact, it has been suggested that non-gaussianity arising during inflation may be useful in probing 
models with a split-SUSY spectrum |146| . In this section we instead focus on the implications of the 
post-inflationary epoch on CMB observations. 


1. Consequences for Inflationary Constraints from the CMB 

CMB data has reached an impressive level of accuracy and can be used to place stringent constraints on 
purposed inflation models. However, to infer which models of inflation are most likely given the data, we 
match observations today to the predictions of the inflation model and this matching requires a knowledge 
of the post-inflationary expansion. The equation of state (e.g. thermal w = 1/3 or non-thermal w = 0) 
during the post-inflationary epoch determines the expansion rate, and thus the rate at which primordial 
perturbations (re)enter the Hubble horizon. For a spectral index n s that is not strictly scale-invariant 
( n s ^ 1) this difference in expansion rate will lead to different predictions for the power spectrum nm- 
1156 1. The matching and the differences that can arise due to two different post-inflationary histories are 
shown in Figure [4] As shown in the figure, given a mode of comoving wavenumber k, we must trace 
back to the instant this mode exited the Hubble horizon (14711155] . This occurs when k = auHk where 
H and a denote the Hubble parameter and scale factor respectively, and a subscript k labels values at 
horizon crossing. Where exp(N(k)) = a en d/ak with N the number of e-folds before the end of inflation 
and a en d the scale factor at the end of inflation. Assuming single-field slow-roll inflation for the purposes 
of illustration and characterizing the post-inflationary expansion by an effective equation of state w. One 
finds the matching equation M 


N(k, w) ~ 71.21-In 
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where p en d is the value of the energy density at the end of inflation, 14 is the inflaton potential as the 
kth mode leaves the horizon, and p r is the energy density at which the universe is assumed to become 
thermalized. The first two terms in (41) are model independent. For GUT scale inflation the third term 
is roughly —10. The fourth term is typically order unity given that the value of the inflaton potential 
necessarily evolves slowly as inflation proceeds. Finally, if the universe thermalizes promptly the last 
term is negligible, and we recover the familiar result that 50 < N < 60 for modes contributing to the 
CMB. 

In [157 


the difference between a SUSY non-thermal history and thermal history was explored and 
using (411 it was shown that the change in equation of state from the thermal case (w = 1/3) to the 


moduli phase leads to a shift 


AN = 


1 — 3w , 

—;-In 

12(1 + w) 



(42) 


with p en( j = 3I4 nc j/2 for a given potential. If w < 1/3 (in fact for the moduli phase we typically 
have w = 0), AN is negative, since p r < p en d with p r the energy density at the end of the moduli 
phase. For simplicity we can assume inflationary (p)reheating was instantaneous and following [ 157] 
we make the substitution p en d —> p a osc , using pf = (7r 2 /30)p*(T/ r )(T/ r ) 4 . At the onset of oscillations 
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FIG. 4: Evolution of physical wavelengths as labelled by their inverse wavenumber k" 1 during inflation 
(below the x-axis) and during the post-inflationary epoch (above the x-axis). The solid (blue) line 
represents the Hubble radius, H^ 1 in a Universe dominated by a radiation fluid w = 1/3, the dashed 
(red) line is the Hubble radius, -ff” 1 in a post-inflationary era dominated by a pressure-less fluid, w = 0. 

We compare the evolution of a physical mode fc* that re-enters at CMB decoupling in the standard 
scenario (Radiation —> Matter —► Dark energy) with a mode k* that re-enters at CMB decoupling in the 
nonthermal scenario (Matter —> Radiation—► Matter —> Dark Energy). These inodes exit the Hubble 
radius at different times during inflation, t* and f 4 , which translates into a shift in the number of e-folds 
AN = HAt. The corresponding shift in the pivot scale or any co-moving mode is given by 

y — h P ~ AN 

Iby. - A/5|cC 
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Posc(t°sc) = ^m 2 Aa 2 and (421 becomes 


AN = —0.04 + ^ In ( 9 *^\ )T ; 

12 V m 2 Aa 2 


= -10.75 + In 
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(43) 


where we used w = 0, and the second line expresses the parameters relative to fiducial values. If scalar 
decay proceeds via a gravitational strength coupling, equation Q eliminates the mass dependence in 
r ^]- With C3 = 1/(47t) we find 


AN = -10.68 + — In 
18 
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(44) 


This shift in the number of e-foldings leads to a shift in both the scalar tilt and tensor-to-scalar ratio 
n „ and r. The uncertainty in is associated with the running a, = dn„/d\nk and to lowest order in 
slow roll [T5gHTB0j 


An s = (n s - 1) 
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The resulting fractional uncertainties A r/r, A?r s /|n s —1| in these observables can be substantial [14811150] . 
In particular, the theoretical uncertainty in n s can be comparable to the precision with which it is 
measured by Planck jlFlj . 

The authors of cm argued that accounting for additional consistencies of models, such as requiring 
a realistic dark matter candidate, could be used to improve the level of constraint. As an example, it 
was shown that working in a SUSY framework, and appropriately accounting for the requirement of 
a successful WIMP candidate could lead to an improved situation for restricting inflation models. As 
we discussed in Section III late decays of the moduli (corresponding to low reheat temperatures which 


maximize AN) will typically dilute thermal WIMPs to an unacceptable level for models with a long 
non-thermal epoch. To remedy this it is usually assumed that WIMPs will be produced non-thermally in 
the decays. However, this implies a large interaction cross-section and large predicted signals for indirect 
detection experiments in cases where all the dark matter is produced in this way. In the MSSM this 
implies Wino or Wino/Higgsino like WIMPs and constraints from FERMI, PAMELA, and in some cases 
direct detection (such as XenonlOO and LUX) can place strong constraints on the cross-section. The 
authors of [1 157] used these constraints to restrict the allowed reheat temperature, which in-turn reduces 
the uncertainties in AN and so the observables in ( |45| ). The bounds on the reheat temperature were later 
improved in Eiuna with the inclusion of data from HESS and taking into account important effects such 
as Sommerfeld enhancement. We will discuss more about non-thermal dark matter and the constraints 
below. Here we want to stress that although allowing for freedom in the post-inflationary history can lead 
to less stringent constraints on inflation model building when only CMB data is taken into consideration, 
properly accounting for other components (and the associated data sets) - in this example we considered 
the origin and properties of dark matter - leads to a new way to begin to probe the so-called dark ages 
prior to BBN. 


2. Dark Radiation and N e ff 

In addition to moduli decays to Standard Model particles and their superpartners, there may also be 
decays to hidden sector fields. Indeed, this is a common expectation in string-based models that give 
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rise to the non-thermal history for dark matter discussed above [162H161] . If the particles resulting from 
decay are light (meaning relativistic) at the time of BBN and/or recombination, and non-interacting 
with MSSM particles, this leads to additional radiation coming from the hidden sector. If these particles 
contribute substantially to the energy density they will affect the expansion rate changing predictions 
for both the abundances of primordial elements EMI and the physics of the CMB EM]- Thus, using 
precision cosmological measurements one can establish constraints on the amount of dark radiation that 
is permitted within a particular class of models - see m and references within. 

During radiation domination (after the decay of the last modulus) the effect of the hidden sector 
radiation on the Hubble expansion can be understood through the Hubble equation ?>H 2 Mp l = p r , where 
the total relativistic contribution to the energy density is 


Pt 3 q 9*T 


(46) 
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where the sums are over relativistic hidden sector particles with degrees of freedom, the factor of 
7/8 results from Fermi-Dirac statistics of fermions, T h is the temperature of the hidden sector particles 
(which importantly need not be equilibrated with Standard Model radiation), and g MSSM is the visible 
sector relativistic degrees of freedom, which in the early universe would be at least g MSSM = 228.75 in the 
MSSM, but near the MeV scale only the photons and neutrinos contribute with 

5mssm(2^) = S 7 T > (48) 


where = 2 for the photon, g„ = 2 for neutrinos, and N e g is the effective number of neutrino species at 
temperature T„. 

At the time of BBN, the neutrino temperature tracks the photons so that T v = T and so with three 
relativistic neutrinos, N e g = 3, the Standard Model prediction is g M ssM = 7.25. However, because the 
neutrinos are weakly interacting (a v v) ~ G 2 F T 2 with Gf ~ 10 -5 GeV ~ 2 Fermi’s constant, they decouple 
below the temperature of BBN (~ MeV) and the entropy in photons increases (so that the total entropy 
is conserved). At the time of recombination we have 

7 

3mssm (Tree) = S7 + gfl'i/ATefif 

7 / 4 \ 4/3 

= 2 + -.2. (3.046) ( n ) 

= 3.385, (49) 



where N e g = 3.046 for the Standard Model m ■ The increase in the photon entropy following neutrino 
decoupling leads to an increase in the temperature of T u /T = (4/ll) 1 / 3 [ 1661 . Given the Standard 
Model (or MSSM) predictions, any new dark radiation would lead to an additional contribution to 3 * 
or N e g. For historic reasons, constraints on new hidden sector radiation are typically expressed through 
iV ef f. In applying these constraints its important to remember that if hidden sector radiation couples 
differently to decaying moduli (or the inflaton during reheating) than Standard Model particles, this can 
lead to different temperatures for each species and this will be preserved in the absence of interactions 
between the systems of particles. As we will discuss soon, such a situation can also lead to isocurvature 
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perturbations. The connection between isocurvature and dark radiation in non-thermal histories was 
studied in |168| . 

Constraints from BBN result from requiring agreement with both the abundances of 4 He and D 
m, which implies N e g = 3.24 ±1.2 at the time of BBN. At the time of recombination, the Planck 
satellite provides constraints with the current results [57| implying N e g = 3.13 ± 0.32 (68% CL), which 
is consistent with A N e ff = 0. 


Dark Radiation and UV Completions. We have seen that axions are a generic prediction of 
string / M-theory and as such light axions can lead to strong constraints in light of bounds on N e f /. 
The consequences have been explored in both the case of M-theory mi Type IIB Large Volume string 
compactifications [211116311170H172] , and more generally |1621116811173H175] . If these axions remain light 
(i.e., if the shift symmetry holds to a good approximation after moduli stabilization), then they could 
constitute dark radiation. Two questions thus naturally arise: (i) how generic is the presence of dark 
radiation in string compactifications?, and ( ii ) if both dark radiation and dark matter have a common 
origin from the decay of a string modulus, can one find a correlation between them by combining dark 
matter observational constraints with present bounds on N e g? These questions have been addressed in 
jl~68l IT681 H76] . 

For question (i), we note that in string compactifications with perturbative moduli stabilisation (like 
the Large Volume Scenarios), the production of dark radiation is unavoidable. This is easy to see, since 
the axionic shift symmetry is protected at perturbative level and broken only by non-perturbative effects. 
Thus, while a modulus gets a large mass during perturbative stabilization, the axionic directions are only 
slightly lifted by non-perturbative effect^] 

With regards to question (ii), in the presence of dark radiation, the reheating temperature T r can 
be written as a function of N e g and the mass of the modulus m a . That is, the modulus decays to 
both visible and hidden degrees of freedom with partial decay widths given by T v ; s = c v i s r 0 and T^m = 
ChidLo, respectively, where c v i s and Chid are model-dependent prefactors and, from Q, we see that To = 
(l/27r)(TO^/TOp). Since the hidden sector is composed of light axions with only gravitational couplings, 
these relativistic decay products do not thermalise and in fact add to the effective neutrino species. 

The expression for A N e ff may be obtained as follows. The conservation of comoving entropy s = 
g(T)a 3 T 3 implies that the temperature and energy density of the visible sector go as T„j S ~ g i/ 3 a an d 

Pvis = i" 3° a 4 ■ Then, at neutrino decoupling, the ratio of the energy densities is given by 
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The number of extra effective neutrino species is then 
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On the other hand, the visible sector reheats to a temperature given by <§, with c replaced by c v - ls . 
Using (51) to eliminate c v i s of ^ in favor of A N e ff and Chid? one obtains the parametric relation between 


19 Light axions can still be removed from the low-energy spectrum by being ‘eaten’ by the anomalous U(l) gauge bosons 
due to the anomaly cancellation. However, this is not the case for light closed string axions living in the bulk. In 
stabilization schemes that dominantly depend on non-perutrbative contributions to the superpotential (like KKLT), the 
axionic direction is strongly lifted. 
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FIG. 5: Constraints on the (AN e g, rncivO-plane for Chid = 1 , g* = 68.5 and m^ = 5 • 10 6 GeV: the solid 
line is based on Fermi data whereas the dashed line represents the freeze-out temperature. The shaded 
region is ruled out due to DM overproduction both in the thermal case (for ?bdm 40 GeV and below 
the dashed line) and in the non-thermal branching scenario (above the solid and dashed lines). 


T r and AN e ff. 


\J AN eff 


(52) 


The constant of proportionality depends on Chid and rn a (for the full expression, we refer to |176| 1. For 
Chid = 1, g* = 68.5 and m 0 = 5 • 10 6 GeV, requiring A N e g < 0.96 translates into a lower bound on the 
reheating temperature: T r h > 0.73 GeV. 


There is one striking consequence of (521. As A N e ff 0, the reheat temperature becomes large 
enough that any dark matter has time to thermalize after being produced. Thus, if future observations 
reveal that the effective number of neutrino species is close to the Standard Model value, the implication 
for dark matter is that a thermal history seems preferred. This is non-trivial information about the 
thermal history of dark matter coming from the CMB. 


One can combine (521 with bounds on the annihilation cross section of dark matter coming from Fermi 
data (we will use the dwarf galaxy constraints described in Section III I. To do this, it is useful to re-express 
the Fermi constraints in terms of the reheat temperature. This can be done as follows. Gamma-ray data 
places an upper bound on the annihilation cross section for given final states as a function of the dark 
matter mass. Assuming that the relic density is satisfied in a non-thermal scenario and is given by 
(101, and using the fact a typical freeze-out temperature is Tf ~ m x o/20 GeV, one can re-express the 
constraints as a lower bound on T r as a function of the dark matter mass. Using (521, this then translates 
into an upper bound on A N e ff as a function of dark matter mass. This can be seen in Figure [5] In 
Figure [6] these same exclusion plots have been displayed for various values of the modulus mass. 
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FIG. 6: Lower bound on the DM mass as a function of AN e g for different values of the modulus mass: 
rricj) = 4 • 10 6 GeV (solid line), ra^ = 5 • 10 6 GeV (dashed line) and m$ = 6 • 10 6 GeV (dotted line). The 
shaded region is ruled out based on Fermi data. Here we have set g * = 68.5 and Chid = 1. 


3. Isocurvature and the CMB 

Consistency with cosmological observations requires that the largest contribution to the temperature 
anisotropy of the CMB must be adiabatic m- Rigid constraints can be placed on models that predict 
departures from strict adiabaticity, however the level of constraint can be model dependent as multiple 
isocurvature components and theoretical priors can play an important role [mj. For a single source of 
isocurvature Planck places strong constraints on model building where isocurvature can be ruled out at 
the percent level m- 

In some cases non-thermal histories with moduli can lead to a departure from adiabaticity and signif¬ 
icant isocurvature contributions. If the fields during inflation have masses comparable to (or less than) 
the Hubble scale, the moduli can fluctuate and lead to a second source of cosmological perturbations - 
differing from the inflaton. This spectrum would be different than that of the inflaton (and its decay 
products) and so depending on the post-inflationary history this could lead to a large source of isocur¬ 
vature perturbations at the time of recombination. For the histories with moduli we have considered in 
this review, an analysis of isocurvature constraints and their importance for non-thermal histories was 
performed in fl68 l (for earlier studies see [ 1781 1T9J). To understand the possible origin of isocurvature in 
non-thermal histories consider the curvature perturbation associated with each type of fluid in the early 
universe (e.g. radiation, dark matter, moduli, etc.). The curvature perturbation of the *th fluid is 

Ci = -* - H ^ . (53) 

Pi 

where 4/ is the Newtonian potential and dots denote time derivatives, and 5pi is the density fluctuation of 
the ith fluid. In connecting with observations it is conventional to define the isocurvature contribution of a 
particular fluid relative to the radiation perturbation so that we define the (gauge invariant) isocurvature 
perturbation of the ith fluid by (e.g. [ 180j ) 


Si — 3 (Ci Cr) ; 


(54) 
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where £ r is the radiation density perturbation. Observables are given by the correlation functions of these 
fluctuations or more conventionally by the power spectrum defined as 

(X(k)Y(k')} = (27r) 3 <5 3 (k - k ')P XY , (55) 

for two fluctuations X and Y. The observables for the isocurvature fraction, on, its correlation with the 
curvature perturbation, r.j, and the cross correlation between any two isocurvature modes, are given 
by 


_ PSjSj 

Ej p SjSj + -Pec ’ 
Ps t C 

y/Ps iSi P« ’ 

= PSjSj 

^/PSiSiPSjSj ’ 


(56) 

(57) 

(58) 


which are all quantities that are constrained by observations m- 

Wh.en do isocurvature perturbations arise in non-thermal histories? If the mass of the modulus is 
lighter than the Hubble scale m a < Hj during inflation, then the quasi-deSitter period will result in long- 
wavelength fluctuations of the modulus field with an average amplitude 6a ~ Hi/2n. This leads to an 
additional source of cosmological perturbations with the corresponding curvature perturbation different 
from that of the inflaton. Following inflationary reheating - where in a thermal history all of the energy 
and matter of the universe is assumed to be created - the modulus can decay leading to an additional 
source of radiation and matter. Thus, whereas radiation and matter created during inflationary reheating 
will inherit the inflaton’s fluctuation £/, those produced from moduli decay will instead be set by the 
modulus (6a) and so (56)-(581 may be non-zero leading to constraints. 

However, even though we have seen that moduli decay is an essential part of a non-thermal history, 
there are still many ways in which isocurvature perturbations may be observationally irrelevant and lead 
to no new constraints on model building. Firstly, if the mass of the modulus is above the Hubble scale 
during inflation m a > Hi, which may be argued to be generic ( 181| 1182] . then in a single Hubble time the 
amplitude of its fluctuations will be exponentially suppressed on large scales by a factor exp(— m 2 a / (3i7|)) 
and so the inflaton will be the only relevant source of cosmological fluctuations (183] . Another important 
observation was made by Weinberg, who demonstrated that even if an isocurvature mode is generated 
initially, if local thermal equilibrium is reached these modes will become adiabatic |184j . And finally, if 
the modulus comes to dominate the energy density of the universe (determining the cosmic expansion 
rate) this also has the effect of washing out any existing isocurvature perturbations. In most cases studied 
thus far arising in string and M-theory negligible isocurvature results 


VI. CHALLENGES AND OUTLOOK 

We have seen that the presence of moduli in the early universe, along with SUSY breaking near the 
TeV scale, motivate a non-thermal history. We discussed how this approach can account for cosmological 
dark matter, axions as a solution to the strong CP problem, and the baryon asymmetry of the universe. 
However, we have also seen that in many cases existing data already puts very tight constraints on each 
of these. The most severe case is for standard SUSY WIMPs EH HH suggesting that more creative 
dark matter model building in the hidden sector may be necessary (see e.g. | l85] ). We have also seen 
that a number of cosmological constraints - particularly those associated with the CMB and structure 
formation - can lead to rigid constraints on these scenarios. As data continues to improve at a rapid pace, 

















33 


this should help us continue to work toward unmasking the cosmological history of the post-inflationary 
universe and shape the direction for model building. 

In most of this review we have emphasized the model independent aspects of non-thermal histories 
arising from moduli. This is possible because the prediction relies on the physics of the hidden sector. 
However, a particular challenge for string/M-theory approaches is embedding the visible sector. This is a 
challenging endeavor, requiring a model dependent approach, but we emphasize here that in cases where 
this is possible and quasi-realistic models can be constructed that the physics is largely independent of 
the hidden sector - a non-thermal history remains. Other model dependent challenges do arise. The 
‘gravitino problem’ is a common example. Gravitinos produced from the decaying modulus can lead 
to both dark matter and dark radiation production. This will depend on model dependent features, 
such as the branching ratio and the gravitino mass. Whether this is a problem or not, depends on the 
model. For example, in the G2 — MSSM, the gravitino mass is comparable to the moduli mass and 
so the decay is kinetically forbidden m- Another challenging problem is the so-called ‘over-shoot’ or 
‘Kallosh-Linde’ problem [45 . As an example, in Type IIB (KKLT) string compactifications, the height 
of the stabilization barrier is set by the gravitino mass. Thus, if the Hubble scale during inflation exceeds 
this barrier height the moduli can overshoot their minimum and de-compactify the extra dimensions. 
This problem is again model dependent and depends on the stabilization mechanism. For example, in 
models of ‘Strong Moduli Stabilization’ the problem can be avoided (120 ;, 1186 , 187 ], A more model 
independent solution may result by having a dynamically changing gravitino mass as discussed in |188[ . 
Model dependent obstructions like the overshoot and gravitino problems emphasize the need for further 
theoretical investigation. Regardless, we have seen that a non-thermal history is a robust prediction of 
hidden sector physics in string / M-theory approaches to beyond the Standard Model physics. 
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